Transport of polar molecules by an alternating gradient guide 
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An alternating gradient electric guide provides a way to transport a wide variety of polar molecules, 
including those in high-field seeking states. We investigate the motion of polar molecules in such 
a guide by measuring the transmission of CaF molecules in their high-field seeking ground state, 
with the guide operating at a variety of switching frequencies and voltages. We model the guide 
using analytical and numerical techniques and compare the predictions of these models to the 
experimental results and to one another. The analytical results are approximate, but provide simple 
and useful estimates for the maximum phase-space acceptance of the guide and for the switching 
frequency required. The numerical methods provide more accurate results over the full range of 
switching frequencies. Our investigation shows that, even when the fields are static, some high-field 
seeking molecules are able to pass through the guide on metastable trajectories. We show that the 
maximum possible transmission requires accurate alignment within the guide and between the guide 
and detector. 

PACS numbers: 37.10.Pq, 37.20.+j, 41.85.Ja 



I. INTRODUCTION 

Cold polar molecules are useful for a wide range of 
applications in physics and chemistry. These include 
tests of fundamental physics [U [2] , high resolution spec- 
troscopy and lifetime measurements [21 U [Sj i controlling 
chemical reactions |B] E] , and implementing quantum in- 
formation protocols El HO] ■ For all these experiments 
it is important to have intense sources of cold molecules, 
and to be able to control their motion. Much progress 
has been made in both directions over the last few years. 
The supersonic expansion method has long been used 
as a versatile technique for producing high-flux beams 
of cold molecules. These beams have been decelerated 
to low speed using the Stark deceleration technique [TT] 
and, more recently, its optical [T? and magnetic [13l [M] 
analogues. The slow beams produced this way have been 
trapped electrically [T5j [TB] and magnetically [T71 [15] . 
A new method for producing very high fluxes of slow- 
moving cold molecules is the buffer gas technique [I9tl20]. 
The molecules of interest are formed inside a cryogenic 
buffer gas cell where they cool to low temperature by col- 
lisions with cold helium. The molecules escaping through 
a hole in the cell can be guided to a room tempera- 
ture, ultra-high vacuum region where the experiments 
take place. The guide separates the molecules from the 
helium and the cryogenic environment while maintain- 
ing the very high density available at the source. Both 
magnetic [5D] and electric guides [H] for low-field seekers 
have already been used for this purpose. In these experi- 
ments the weak-field seeking molecules were confined by 
field minima created on the axis of the beam line. Since 
a local field maximum cannot be created in free space 
the same methods cannot be used to guide a beam of 
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high- field seeking molecules j^^. There is a strong 
motivation to develop guides for high-field seekers, since 
all ground-state molecules are high-field seekers and all 
heavy molecules in low-lying quantum states seek high 
electric fields for all practically relevant field strengths. 

To guide molecules in high-field seeking states the 
method of alternating gradient (AG) focussing can be 
employed. The molecules are guided as they pass through 
a sequence of lenses, each of which focusses in one trans- 
verse direction and defocusses in the other, the direc- 
tions alternating from one lens to the next. They are 
guided because, in each direction, they tend to be fur- 
ther from the axis in the focussing lenses than in the 
defocussing lenses, and the applied forces are approxi- 
mately linear in the off-axis displacement. This tech- 
nique can be used to selectively guide molecules in either 
high-field or low-field seeking states. It was first applied 
to charged particles p4l and later extended to neutral 
molecules [25, 26, 27 . More recently, prototype AG de- 
celerators have been demonstrated [251 US ISHl 131] and an 
AG guide has been used to separate the two conformers 
of CgHyNO by using the guide's sensitivity to the ratio 
of dipole moment to mass |32]. An AG structure with 
a double-bend was also used to guide and filter the slow 
fraction of an effusive beam [33], but the high- and low- 
field seekers could not be distinguished in these experi- 
ments since the detection method was not state-selective. 
Focussing of large molecules, with particular emphasis on 
state- and conformer-selection, is discussed in [34] . 

In this paper, we demonstrate AG guiding of a super- 
sonic beam of CaF molecules. Since our detector resolves 
the rotational and hyperfine structure we can investigate 
the effect of the guide on individual quantum states, and 
we choose to focus on the high-field seeking ground state. 
We measure the efficiency of the guide as a function of 
the applied voltage and the switching frequency. We 
also model the guide using both analytical and numer- 
ical techniques, showing how the molecular trajectories 
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FIG. 1: (Color online) (a) Setup of the experiment, (b) 
The central region of the guide showing the contours of elec- 
tric field magnitude in each of the two switch states. The 
magnitude of the electric field, in kV/cm, is given by the 
color/grayscale map for the case where V=10kV. The cen- 
tral electrodes of the guide are shown in grey and are labelled 
A,B,C,D. Their polarities in each of the two switch states are 
indicated. 

and the phase-space acceptance of the guide change as we 
move from an ideahzed model where the forces are linear 
functions of position, to one that includes the most im- 
portant nonlinear forces, and then finally to a complete 
numerical simulation. We compare our experimental and 
theoretical results. Finally, we demonstrate that some 
high-field seeking molecules are focussed to the detector 
even when the guide fields are static. This surprising 
result is due to molecules that are first ejected into the 
outer regions of the guide but later pulled back in towards 
the centre. 



II. SETUP OF THE EXPERIMENT 

The experimental setup is shown in figure [ija) . The 
methods we use to produce and detect cold supersonic 
beams of CaF are described in [5S] , so we give only some 
brief details here. A pulsed solenoid valve containing a 
4 bar mixture of 98% Ar and 2% SFg emits gas pulses 
of « 200^s duration into a vacuum chamber. When 
operating at a repetition rate of 10 Hz, the pressure in 
this source chamber is typically 10^'^mbar. A Ca tar- 
get held outside the valve nozzle is ablated by a pulsed 
Nd:YAG laser of wavelength 1064 nm, energy 25 mJ, and 



10 ns pulse duration. The hot cloud of ablated Ca reacts 
with the SFg to form CaF radicals which are then en- 
trained in the expanding carrier gas and cooled to «3 K. 
In our coordinate system, the valve nozzle is at the ori- 
gin, the z-axis is the guide axis, the x-axis is parallel to 
the line connecting the centres of electrodes A and B, and 
the ablation laser fires at i = 0. The CaF beam passes 
into a differentially pumped chamber maintained below 
lO^^mbar, via a 2 mm diameter skimmer at z = 70 mm. 
The molecules pass through the guide whose front end is 
at z = 125 mm, then through a 5.8 mm diameter aperture 
ai z — 1160 mm, and are finally detected at z — 1215 mm 
by time-resolved laser-induced fiuorescence (LIF) using 
a cw dye laser at 606.3 nm. The 220 iiW probe laser 
beam has an approximately uniform intensity distribu- 
tion within a 2.8 mm diameter circle, crosses the molecu- 
lar beam at right angles and makes an angle of 45° to the 
X-axis. The ablation laser beam is parallel to the probe 
laser beam. 

The guide is formed by four parallel cylindrical stain- 
less steel rods of length Im and diameter 6 mm, with 
spherically-rounded ends of radius 3 mm. These rods are 
arranged with their centres on the corners of a square of 
side 7.87 mm, and their axes parallel to the beamline. To 
make sure these rods are supported rigidly, each one is 
attached to one of four larger rods, Im long and 16 mm 
in diameter, using metal dowels at five equally-spaced po- 
sitions along their lengths. These larger rods have their 
centres on the corners of a square of side 33.9 mm. To 
guide high-field seeking molecules the guide is switched 
between the two states shown in figure [T|[b), by reversing 
the polarity of electrodes B and C using a bipolar high 
voltage switch. The voltages applied to electrodes A and 
D are fixed. In switching between the two states the 
saddle-shaped field is rotated through 90° . In state 1 the 
molecules are focussed along x and defocussed along y, 
while in state 2 the focussing and defocussing directions 
are reversed. 

We determine the number of ground-state CaF 
molecules transmitted to the detection region as a func- 
tion of the two main experimental parameters, the ap- 
plied voltage, V, and the switching period, 2r (see Fig.jlJ. 
To obtain each dataset we fix the value of these guiding 
parameters and then scan the frequency of the probe laser 
over the two hyperfine components of the X^S+(ti" — 
0) - A^Ili/2{v' = 0) Q(0) line. Then, we step one of 
the guiding parameters and repeat. To remove the effect 
of the source flux varying from one dataset to the next, 
we modulate, on a shot-by-shot basis, between two val- 
ues of the guiding parameters, one of which is held flxed 
throughout and used as a constant reference to normal- 
ize all the other data. The amplitude of the two spectral 
lines determines the number of ground-state molecules 
passing through the detector |35j, and their widths are 
related to the angular divergence of the beam exiting 
from the guide. 
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III. MODELLING THE GUIDE 

The theory of AG guiding was developed in the context 
of charged particle acceleration [21] . The application of 
these ideas to the focussing of neutral particles was first 
given in [2S], and later extended in [301 123 • Here, we 
apply the theory to our four-rod guide and show how 
the results are modified as various approximations are 
removed. 

The development of the theory begins with a suitable 
multipole expansion of the electrostatic potential of the 
guide ^{x,y), as a function of the transverse coordinates 
X and y: 



^{x,y) = $0 Ui- 
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Here, tq characterizes the transverse size of the guide, <i>o 
the applied voltage, and oi, 03, 05 the relative sizes of 
the dipole, hexapole and decapole terms in the expansion. 
Assuming that 05 ^ 03 <g; oi the electric field magnitude 
in the region where x,y <rQ is well approximated by 

E{X, Y) = Eo{l + 03(^2 - y2) _^ (2a^ - 6a5)X'^Y^ 

+ a,{X^ + ¥') + ■■■), (2) 

where Eq is the electric field magnitude on the axis of 
the guide. We have set oi = 1, introduced the scaled 
coordinates X = x/ro, Y = y/ro, and neglected terms of 
order 0305, 03 and higher. We expand the Stark shift in 
a Taylor series about Eq, 

W{E) = W{Eo) ~ tiMEo)iE -Eo) + -- - , (3) 

where fJ,cs{Eo) — —^\eo is an effective electric dipole 
moment whose value converges to the molecular dipole 
moment in the limit of strong rotational mixing. We ne- 
glect the higher-order terms in the Taylor expansion be- 
cause, in the strong fields of the guide, the Stark shift of 
the molecule is nearly linear in the electric field strength. 
The force on the molecules is the spatial gradient of their 
Stark shift. Thus we arrive at the following approximate 
equations of motion: 

x"^n^(^x + xY^x + cxA, 



Y" = 



0-3 
l«3| 



Y + xX^Y + (Y' 



For a four-rod guide such as the one demonstrated 
here, the ratio of rod radius to rod separation has a large 
effect on the efficiency of the guide. When the rods are 
large the forces are weak and the confinement is poor, but 
if the rods are made too small the nonlinear terms in the 
force become too large and the confinement is again poor. 
These details were examined for a general AG guide in 
|36j and the best values of 03 and 05 were determined. To 
obtain values of Eq, 03 and 05 that describe our guide, 
we calculated the electric field in the central region of the 
guide using a finite element model and then fitted this to 
Eq. Q. The electric field is shown in Fig.[T](b) and the 
best fit values are Eq = 4.395ykV/cm, 03 = 0.120 and 
05 — 0.0022, where V is the applied voltage in kV. From 
these, we obtain ( = 0.0363 and x = 0.131. We note that 
the values of 03 and ( are close to those recommended in 
[36] for optimizing the guiding efficiency (see in particu- 
lar Fig. 12 of [36]). The value of depends on the applied 
voltage and the molecular mass and dipole moment. In 
this section we study ground state CaF molecules trav- 
elling through our guide with the applied voltage set to 
8kV, which gives H. = 7776rads~^. 

Let us now calculate the molecular trajectories and the 
phase-space acceptance of the AG guide. We find it useful 
to distinguish three levels of calculation. At the simplest 
level, the nonlinear terms in the equations of motion are 
neglected so that the trajectories are given by ^ with 
X = C = 0. The solid line in Fig.|2ja) shows an exam- 
ple trajectory obtained within this linear approximation 
for the case where the switch duration is r = 160 /is (i.e. 
Q.T ~ 1.24). The trajectory is characteristic of AG fo- 
cussing. It consists of an oscillation at the driving fre- 
quency with a relatively small amplitude - the micro- 
motion - superimposed on a lower frequency oscillation 
of larger amplitude - the macromotion. In this linear 
model, molecules are only transmitted by the guide if 
their trajectories fit inside a square aperture of side 2xa 
whose size is chosen to reflect the natural boundaries set 
by the guide's electrodes. For our guide, we take this to 
be the square whose corners touch the four electrodes; it 
has xa = 1.8mm and is shown in Fig.[ljb). The ana- 
lytical result for the acceptance of the guide is given by 
equation (17) of [36]^: 



^linear — 



7r2rj22,4 [j^ _ cos2(rjT) C0sh^(r2T) 

[cosh(riT) sin(rir) + sinh(rir)]2 



(6) 



(4) 



The solid line in Fig.[2];b) is a plot of this formula as a 
function of r, when fl — 7776rads~^. When r is very 
short the micromotion amplitude is very small and the 



Here, differentiation is with respect to time, ri/(27r) is 
the oscillation frequency in a focussing lens, and x and 
C represent the strength of nonlinear coupling and cubic 
terms in the force. These parameters are given by 
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^ In this paper, we use the acceptance in (x,x',y,y') space. In 
cases where the transverse coupling is neglected, this can be 
written as the square of the acceptance in just one direction, 
say {x,x'), which is how this result was written in |36| . Here we 
have also included the size of the aperture xa, which was set to 
ro in |36l . 
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FIG. 2: (Color online) AG guiding of ground state CaF 
molecules calculated using various levels of approximation, 
(a) An example trajectory for a molecule with initial condi- 
tion x{0) = 1/(0) = 0, v^{0) = 1.8m/s, Vy{0) = -1.2m/s, 
Vz — 600 m/s. The guide parameters are V = 8kV and 
r = 160 us. (b) Phase-space acceptance as a function of r. 
The solid red lines are calculations based on a linear force 
with Q = 7776 rads^^. The dashed green lines include the 
cubic term in the force, with coefficient ( — 0.0363. The dot- 
ted blue lines (connecting blue circles in (b)) include both the 
cubic term and the coupling term, with coefficients ^ — 0.0363 
and X ~ 0.131. The dot-dashed purple lines (connecting pur- 
ple squares in (b)) are obtained from a complete numerical 
simulation using the exact electric field map, Stark shift and 
electrode geometry. 



converging effect of each focussing lens is almost exactly 
undone by the divergence in the following defocussing 
lens. As a result, the macromotion wavelength is very 
long and the acceptance is small. As r increases so too 
does the micromotion amplitude. The molecules tend 
to be considerably closer to the axis in the defocussing 
lenses than in the focussing lenses, and so the latter have 
a greater effect resulting in an increased acceptance. In 
our example, the acceptance has a maximum value of 
350(mmm/s)^ when r ~ 160 /iS. If r is increased fur- 
ther the acceptance rapidly decreases as the molecules 
are overfocussed. The motion is unstable for t > 240 /is 
and the acceptance falls to zero. 

The next level of calculation uses the full equations of 
motion, Eq. to calculate the trajectories of a set of 



molecules numerically. The dashed line in Fig.|2ja) shows 
how the trajectory is modified when the cubic term is 
added to the force. Since ( is positive, this weakens the 
focussing lenses and strengthens the defocussing lenses, 
causing the molecules to deviate further from the axis 
than before. This naturally decreases the phase-space 
acceptance. Just as in the linear case, a molecule is only 
transmitted if its trajectory fits inside the aperture whose 
size is chosen to reflect the natural boundaries of the 
guide. Provided the guide is not operated too close to the 
stability boundary (r — 240 /is in this case) there is an 
approximate analytical result for the phase-space accep- 
tance that includes the effect of the cubic term in ([4| but 
neglects the coupling term. This result, which is equa- 
tion (42) of Pl], is shown by the dashed line in Fig.[2j^b). 
As expected, we see that the cubic term reduces the ac- 
ceptance. The reduction is greatest for small lens lengths 
where the net focussing is already weak and the further 
weakening caused by the nonlinear term has a great ef- 
fect. Because of this, the peak of the acceptance curve is 
shifted to longer lens durations, and its value reduced to 
200 (mmm/s)^. Next, we add in the effect of the nonlin- 
ear coupling term with coefficient x = 0.131. This force 
both weakens the conflnement of the molecules along the 
diagonals of the guide (x is positive) and couples together 
the motion in the x and y directions. The modified tra- 
jectory is given by the dotted line in Fig.[2|^a), showing 
that the macromotion amplitude has increased even fur- 
ther. The phase-space acceptance now has to be calcu- 
lated numerically, and the results obtained are shown by 
the circular points in Fig.|2jb). For small r, and values 
of r close to the cutoff, the acceptance seems to have in- 
creased. This is because the numerical calculation is for a 
Im long guide, while the approximate analytical results 
given in [35] apply to infinitely long guides. Because the 
guide has a finite length some molecules are able to pass 
through even when their trajectories would ultimately 
not be transmitted. At intermediate lens durations the 
acceptance is reduced by the coupling term as a result of 
the overall weakening of the confinement. The optimum 
value of T is again shifted to larger values and the peak 
acceptance shghtly reduced. 

The nonlinear terms play a large role in limiting the 
phase space acceptance and so their inclusion provides 
a far more accurate determination of the guide's accep- 
tance than obtained in the linear approximation. This 
description is still incomplete however because Eq. (|4| is 
only accurate close to the axis, the Stark shift is not 
perfectly linear in \E\, and the real boundaries of the 
guide are set by the electrodes rather than by an artifi- 
cial aperture. The third level of calculation involves the 
complete numerical modelling of the experiment, using 
the electric field obtained from a finite element model, 
the correct dependence of the Stark shift on \E\, and the 
actual boundaries set by the electrodes. The dot-dashed 
line in Fig.|2ja) shows the trajectory calculated using this 
complete numerical model, with the same initial condi- 
tions and guide voltage as before. This trajectory follows 
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closely the one obtained from Eq. Q , demonstrating that 
the simple application of Eq. Q is quite accurate in this 
case. The full simulation predicts slightly larger macro- 
motion amplitudes because the higher-order terms not 
included in Eq. Q tend to further reduce the confine- 
ment, particularly at large distances from the axis. We 
note that if we choose a far smaller initial divergence so 
that the molecule remains close to the axis, the trajecto- 
ries calculated using the different methods all converge. 
The square points in Fig.|2]^b) show the acceptance cal- 
culated using the full numerical simulations. In these 
simulations, we included the final aperture placed down- 
stream of the guide in the experiment. This is important 
because our simulations show that even for lens dura- 
tions well beyond the cutoff, a large number of molecules 
can reach the end of the guide on trajectories that go 
very far from the axis but do not crash into any elec- 
trodes. We study this in more detail in section |V] Only a 
small fraction of these molecules are able to pass through 
the final aperture, and so when included we recover the 
cutoff near the expected lens duration. Over all, the ac- 
ceptance obtained from these full numerical simulations 
closely follows the results obtained using Eq. ^, again 
showing that this simpler (and faster) method of calcu- 
lation provides rather accurate results. 

Figure ^h) illustrates the importance of choosing the 
optimum switching period for the guide. The nonlinear 
terms shift the optimum period upwards to a value very 
close to the cutoff length predicted in the linear theory. 
When the period is too short the nonlinear terms greatly 
reduce the acceptance, but the peak acceptance is only 
reduced by a factor of 2 from that predicted by the linear 
theory. This observation is important because it shows 
that a well designed guide can offer an acceptance that 
is similar to that promised by simple calculations. 



IV. AG GUIDING RESULTS 

Figure [3] shows the ground state CaF signal detected in 
the experiment as a function of r for three different values 
of the applied voltage, 9.5 kV, 8kV and 5.5 kV. In each 
case the signal is the amplitude of the F = 1 hyperfine 
component in the laser-induced fluorescence spectrum, 
and these amplitudes are normalized to the equivalent 
amplitude observed when the guide is turned off. To min- 
imize the effect of source flux variation, the modulation 
procedure described in Sec.[ll] was used for the 5.5 and 
9.5 kV datasets. It was not used for the 8kV data which 
explains why the error bars are larger here. For all the 
data, the central speed and the longitudinal translational 
temperature were determined from the distribution of ar- 
rival times and found to be 600 m/s and 3 K respectively. 
For reference = 8640,7776 and 6120 rads"^ for the 
three voltages used here. In the same flgure we show the 
results obtained from full numerical simulations of the 
experiment. We first discuss the experimental data and 
then turn to the comparison with the simulations. 
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FIG. 3: (Color online) Ratio of fluorescence signal with guide 
on and off versus lens duration for three applied voltages: 
(a) 9.5 kV, (b) 8.0 kV and (c) 5.5 kV. Lines show the corre- 
sponding results obtained from numerical simulations of the 
experiments. Dashed lines are for perfect alignment of the 
probe laser beam and solid lines are for a 2 mm misalign- 
ment. Note the different scales used on the vertical axes of 
the three graphs. 

As expected, the number of guided molecules is small 
when the lenses are short. In the idealized picture of an 
AG guide, in which the focussing and defocussing powers 
are equal, the guide does nothing at all when the switch- 
ing period is very short because the defocussing impulse 
in one lens is almost exactly undone by the focussing im- 
pulse in the next lens, and so the motion-averaged force is 
very close to zero. In the real AG guide, when the period 
is short the nonlinear contributions to the force result in 
a net defocussing of molecules that deviate too far from 
the axis. These molecules are forced out of the guide and 
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the number of molecules reaching the detector is smaller 
than obtained when the guide is off. As the switching 
period is increased the alternating gradient guiding be- 
comes effective and so the number of molecules reaching 
the detector increases. If the lenses become too long, the 
motion becomes unstable as the molecules are overfo- 
cussed and the number of transmitted molecules rapidly 
falls to zero. Increasing the applied voltage shifts the op- 
timum switching period to shorter values such that the 
corresponding value of f2r remains approximately con- 
stant at about 1.5. Increasing the applied voltage from 
5.5 kV to 8kV increases the observed signal by a factor 
of 1.6, which agrees very well with the H.^ scaling ob- 
tained in the linear theory (see Eq. ([6])). However, this 
same scaling predicts that increasing the voltage from 
8kV to 9.5 kV should increase the acceptance by a factor 
of 1.2, whereas in our data the peak signal at 9.5 kV is 
approximately the same as that at 8 kV. 

To simulate the experiments we calculate the electric 
field of the guide using a finite element model and use this 
to determine the position dependence of the molecule's 
Stark shift. The applied force is the spatial gradient of 
this Stark shift. Apart from the regions near the two 
ends of the guide, the electric field is not a function of z. 
In simulating the trajectories, we can choose to neglect 
the field variation at the two ends so that the simulation 
is simpler and faster, or we can include it so that the 
results are more accurate but are slower to obtain. We 
have done both, and find that the proper inclusion of the 
ends has little effect on the results obtained. They are 
included in all the simulations presented in this section. 
Molecules are propagated from the source to the entrance 
of the guide via the skimmer, their trajectories through 
the guide are calculated, and they are finally propagated 
through the downstream aperture to the plane of the de- 
tector. For each value of V and r we determine the set 
of molecules that reach the detector without hitting any 
surfaces along the way. To complete the simulation we 
convert the molecule number to a photon number by cal- 
culating the detection efficiency for each molecule, which 
is a function of its position, its velocity component along 
the probe laser and its transit time through the laser, us- 
ing a simple three-level rate model along with the known 
laser intensity distribution, branching ratio and natural 
linewidth |35l I37j . Since the distribution of transverse 
positions and velocities is different for the guided and 
un-guided beams, the detector efficiency does not simply 
divide out in the signal ratio between guide on and off. 

The results of these simulations, expressed as the signal 
ratio between guide on and off, are shown by the dashed 
lines plotted in Fig. [3] The simulations accurately pre- 
dict the optimum value of r for all three voltages, and 
have the same qualitative shape as the experimental data, 
but predict a considerably higher signal ratio than is ob- 
served experimentally. There are a few possible reasons 
for this discrepancy. Firstly we note that we have only 
partial information about the source. The time of flight 
profile gives the distribution of forward velocities, but the 



distribution of transverse velocities and positions in the 
source is unknown. In the simulations presented here, 
we have assumed normal distributions in all dimensions, 
with a full width at half maximum of 4.7 mm in x and 
y and of 38m/s in Vx and Vy. Decreasing these widths 
tends to decrease the signal ratios predicted by the simu- 
lations, reducing the discrepancy with the experimental 
results, because a smaller, less divergent source increases 
the relative number that reach the detector when the 
guide is off. We think it unlikely that this could account 
for the full discrepancy since the source would need to 
be unfeasibly small or unusually well collimated. Mis- 
alignments of the guide's electrodes also tend to reduce 
the transmission of the guide and so could also be partly 
responsible for the observed discrepancy. Misalignments 
were thought to be the cause of similar discrepancies of 
signal amplitude in experiments with AG decelerators 
e.g. [28]. Finally, misalignment of the probe laser from 
the z-axis dramatically reduces the signal ratio because 
the distribution of molecules at the detector tends to be 
considerably smaller when the guide is on than when it is 
off. The solid lines in Fig.|3]shows the results of the same 
simulations when a probe misalignment of 2 mm is intro- 
duced. This misalignment brings the simulation results 
into far better agreement with the experimental results 
for all three data sets. In the experiments, the probe 
beam could not be moved far in the y direction without 
increasing the background laser scatter to an unaccept- 
able level. This constraint on the probe position makes 
a misalignment quite feasible and we think it is the most 
likely cause of the observed discrepancy. This hypothesis 
is further corroborated in the next section. 

On the high-T side, where the stability cutoff occurs, 
the simulations predict that the cutoff should be at 
slightly higher values of r than observed in the exper- 
iments. This is true of both the simulations with the 
misaligned probe laser and the simulations with perfect 
alignment scaled down so that the peak signal matches 
the data. When the guide is operated close to this stabil- 
ity cutoff, it becomes particularly sensitive to misalign- 
ments of the electrodes since these tend to destabilize 
trajectories that would otherwise be stable. To see this 
effect, we have included one type of imperfect alignment 
in our simulations. We suppose that the guide is perfectly 
straight and that the geometry is perfectly homogeneous 
along the length except for an overall scaling of the en- 
tire geometry which is allowed to vary along the length. 
This particular variety of imperfection, though unrealis- 
tic, has the virtue of being easy to simulate. To quantify 
the scaling required to best follow the misalignments of 
our guide, we measured the size of the gaps between elec- 
trodes A & B, B & D, D & C and C & A (see Fig.[l]), at 
many different positions along the length, and used this 
information to assemble a ^-dependent geometry-scaling 
function. The variation in the gap size arises from the 
stresses imposed at the five support positions. The small- 
est gap size is at the centre, where it is approximately 
5% smaller than the mean value, while at the two ends. 
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FIG. 4: (Color online) Normalized signal versus DC applied 
voltage. The dashed line shows the corresponding simulation 
results for perfect alignment, while the solid line has the probe 
beam misaligned by 2 mm. 

the rods, being unsupported, tend to splay outwards, in- 
creasing the gap size by about 15% over the mean. The 
simulation results presented in Fig.|3]already include this 
imperfection. If we do not include it, the discrepancy 
between experiment and simulation regarding the cutoff 
position increases, and the optimum value of t shifts to 
higher values. The acceptance calculation shown by the 
dot-dashed line in Fig.[2][b) did not include either the 
guide ends or the variation in gap size, which explains 
why the optimum r seen in that figure is higher than in 
Fig. [3j These observations, along with the highly sim- 
plified nature of the imperfection that we have included, 
leads us to speculate that the remaining discrepancy ob- 
served in Fig. [3] is indeed the result of electrode misalign- 
ments. 



V. DYNAMICS OF HIGH-FIELD SEEKING 
MOLECULES IN STATIC FIELDS 

In this section we show that some ground state CaF 
molecules are able to pass through the guide even when 
the fields are static. Guiding of strong-field seeking 
molecules by static fields has been demonstrated previ- 
ously, either by sending the molecules on circular orbits 
around a central wire |38j . or by placing on the beam axis 
electrodes that are so small that the molecules tend to 
miss them [33] • In these cases, the axis of the guide is not 
charge-free, and a static maximum exists there. In the 
present experiments the situation is quite different since 
there are no electrodes on the axis of the guide, and the 
field on the axis is a saddle-point. 

To investigate static guiding, we measured the signal 
as a function of V without switching the guide at all. In 
this case the guide would seem to be a a single 1 m-long 
lens, focusing in one transverse direction, and defocusing 
in the other. Figure |4] shows the normalized ground state 
CaF signal obtained as a function of the applied DC volt- 
age when the speed is 600 m/s. The linear analysis of a 



very long AG guide suggests that no molecules will be 
transmitted once the voltage exceeds V = 0.74 kV, but 
this is not what happens. As the voltage is increased 
from zero the signal ratio does drop, but not as rapidly 
as the linear theory would suggest. Instead of dropping 
to zero it reaches a minimum value of 0.23 at 2.5 kV. The 
signal then increases with voltage, reaching a maximum 
value of 0.64 at 5.5 kV before dropping again. A similar 
surprising behaviour is predicted by the simulations. The 
dashed line in Fig.|4]shows the results of simulations done 
in the same way as described in Sec. |IV[ except that the 
end-fields and the position dependent gap size were not 
included. The distribution of molecules at the source was 
the same as before. The simulations agree well with the 
data between and 2.5 kV, but not at higher voltages. 
The greatest discrepancy is between 3 and 5 kV and we 
find that at these intermediate voltages the transverse 
spread of the molecules is small in the plane of the de- 
tector, and so the results become sensitive to the exact 
position of the probe laser. We find far better agreement 
between simulation and experiment when, in the sim- 
ulation, the probe laser is displaced by 2 mm from the 
z-axis, as shown by the solid line in Fig. |4j This is the 
same misalignment suggested by the AG data discussed 
above. 

Figure |5ja) shows how high-field seeking molecules 
reach the end of the guide when the fields are static. The 
figure shows a map of the electric field magnitude within 
the central 24 mm x 24 mm region of the guide. We see 
that there are three saddle-shaped regions, the one in the 
centre typically considered for AG guiding, and one be- 
tween each of the two positive-negative electrode pairs. 
The figure shows an example of the motion of a molecule 
that starts at the origin and is travelling in the positive- 
y direction. It is accelerated towards the field maximum 
located between the rods, decelerated again as it over- 
shoots and moves into the outer regions of the guide, 
then attracted back towards the centre repeating its mo- 
tion in the negative-?/ direction. While in the vertical 
direction molecules such as this one are defocussed when 
near the guide centre and focussed when in the outer 
regions, the opposite is true in the horizontal direction. 
As a consequence of their motion, these molecules see an 
alternating focus-defocus sequence even though the fields 
are static. In the horizontal direction the average force 
is far weaker and so the stable trajectories are all elon- 
gated, similar to the one shown in the figure. For the 
molecules to avoid hitting the electrodes, they need to 
have a relatively high initial velocity along y. If this is 
too small, the molecules spend too long in the region be- 
tween the two oppositely charged electrodes and tend to 
crash into those electrodes. This leads to an unusual sit- 
uation in which molecules with large transverse velocities 
are transmitted to the end of the guide, while those with 
smaller transverse velocities are not. We have plotted the 
trajectory followed during the first 3.7 ms which is more 
than double the time it takes for a 600 m/s molecule to 
pass through our Im guide. 
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locities that are mostly smaller than this, meaning that 
the majority of the molecules are unable to escape the 
guide. They are lost because they hit the electrodes, not 
because they leave the guide. Figure [5](c) shows y{z) for 
a range of applied voltages, each with the same initial 
condition. We see again the vertical oscillation shown 
in (a), and we see how the wavelength of this oscilla- 
tion decreases with increasing voltage. At ^ 6kV the 
molecule is back at the centre when it leaves the guide. 
The oscillation wavelength does not change much with 
amplitude and so there is a set of stable molecules that 
are focussed back onto the axis at this particular volt- 
age. These molecules are then able to pass through the 
exit aperture and into the detector. As the voltage is 
increased, the beam is focused earlier in the guide and 
exits with larger divergence so that fewer molecules can 
pass through the aperture to the detector. This explains 
the behaviour observed in Fig.|4j 

When V « 12 kV the molecules execute one complete 
oscillation inside the guide and exit close to the axis 
again, as shown in Fig.js] (c). However, this does not 
lead to a second increase in the fluorescence signal as we 
approach this voltage. Rather, as V increases we see a 
considerable broadening of the spectral line which splits 
into two resolved components when V > 8 kV. This split- 
ting comes about because the transmitted molecules all 
tend to have large values of | Uj, | , as explained above. Con- 
sequently, the molecules exiting the guide have a distri- 
bution of Vy values consisting of two peaks either side 
of zero. The separation of the two peaks grows with in- 
creasing voltage, and this is reflected in the distribution 
of Doppler shifts that we measure. 



FIG. 5: (Color online) (a) A contour and color/grayscale map 
of the electric field magnitude in the guide. The field is in 
kV/cm when V — 10 kV. The grey circles indicate the elec- 
trodes. The supporting rods are just visible in the corners of 
the plot. The bold line shows the trajectory, projected onto 
the a;y-plane, of a molecule that is guided when the DC volt- 
age is 6 kV. In this example, the molecule was launched from 
the centre of the guide with = 1 m/s and Vy = 13 m/s. The 
trajectory followed for the first 3.7 ms is shown, (b) The Stark 
shift of ground state CaF as a function of position along the x 
(solid) and y (dashed) axes of the guide, for V = 8 kV. (c) Tra- 
jectories, y versus z, for applied voltages of 2, 4, 6, 8, 10 and 
12 kV. In each case the molecule was launched from the centre 
of the guide with — 0, Vy — 12.5 m/s and Vz = 600 m/s. 



Figure |5] (b) shows the Stark shift of the molecules 
as a function of position along the principal axes of 
the guide when V — 8kV. A molecule that enters the 
guide at the origin will be pushed out along y into the 
outer regions of the guide, but it cannot escape from the 
guide completely unless its initial transverse kinetic en- 
ergy is greater than /ix26GHz (1.25 K in temperature 
units). For CaF this corresponds to a transverse velocity 
of 19 m/s. The molecules are geometrically constrained 
to enter the guide near the centre and with transverse ve- 



VI. CONCLUSIONS AND DISCUSSION 

In this paper we have thoroughly investigated the AG 
guiding of CaF molecules in a high-field seeking quantum 
state. The methods we use are applicable to a wide range 
of other polar molecules 



III we outlined various 
from an analytical 



In Sec. 

methods of modelling an AG guide, 
theory based on linear forces through to full numerical 
simulations. We showed how the trajectories and phase- 
space acceptance change as we make the model more 
realistic. The various levels of modelling are all useful 
and the method of choice depends on the accuracy re- 
quired. For our guide, the peak acceptance predicted by 
the linear theory is an overestimate by a factor of about 
2 compared to the numerical results - this is useful for 
estimating the peak acceptance, though it is important 
to note that the overestimate is larger for short values of 
T. The analytical result for the acceptance that includes 
the cubic term (dashed green hue in Fig.[2];b)) IS a par- 
ticularly useful result since it can be evaluated quickly 
and gives a result that is in fairly good agreement with 
the acceptance calculated numerically. In Sec. |IV| we 
demonstrated the AG guiding of high-field seeking CaF 
molecules, showing how the guiding efficiency varies with 
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switching period for various values of applied voltage. To 
understand the results we simulated the experiments and 
found that the predicted shape of the transmission curve 
agrees well with the data, except near the long-period 
cutoff. By including a simple type of imperfect alignment 
in these simulations, we find evidence that misalignments 
of the guide are responsible for the discrepancy near the 
cutoff. The absolute magnitude of the signal was not in 
agreement with the data, most likely due to a misalign- 
ment of the detector which we could not easily correct. 
In Sec. |V] we demonstrated that a substantial number 
of molecules are transmitted through the guide to the 
detector even when the fields are static, because some 
molecules are first pulled out into the outer regions of 
the guide and then accelerated back towards the centre. 

In presenting our data we have normalized all signals 
to the signal obtained with the guide turned off. This 
ratio provides a useful comparison, but it is important 
to understand that it will depend strongly on the tran- 
sit time of the molecules through the guide, T. When 
the AG guide is on, the number of molecules transmit- 
ted should be independent of T, provided this is long 
enough for the molecules to complete one macromotion 
cycle. But when the guide is off a molecule is only trans- 
mitted if its transverse velocity allows it to pass through 
the guide's aperture, and for a fixed transverse velocity 
distribution the number of these molecules falls as l/T^. 
It follows that the ratio we measure should increase if T 
increases i.e if the guide is lengthened or the forward ve- 
locity is reduced - guides are most useful for transporting 
slow molecules over long distances. The AG guide used 
here could be very usefully employed in conjunction with 
a cryogenic buffer gas source so that molecules in high- 
field seeking states can be guided away from the buffer 
gas to a region of high vacuum. In this case, the molecules 
will enter the guide at considerably lower speeds than in 
the experiments described here. To make a first estimate 
of how well this might work, we consider the use of our 
guide to transport CaF molecules away from a 4 K buffer 
gas cell having a 1 mm diameter opening. We suppose 
that the beam exiting the cell is approximately effusive 
but, as is common [T2j, has a forward velocity-boost that 
shifts the peak velocity to ^ 80m/s. We use a Gaussian 
distribution of transverse velocities with a full width at 
half maximum of 47m/s, place the entrance of the guide 



8 mm from the source, place a 1 mm diameter aperture 
5 mm downstream of the guide's exit, and set V = 8kV 
and T = 200 /iS. A full numerical simulation of this setup 
that includes the z-dependence of the field at the ends of 
the guide shows that of all the molecules emerging from 
the cell, about 0.5% will be guided through the final aper- 
ture. Laser ablation yields up to 5 x 10^^ ground state 
CaF molecules per pulse HO]. With the buffer gas cell 
operated in the partially hydrodynamic regime, we can 
assume that at least 5% of all these molecules enter the 
beam (40% has been demonstrated [ID])) and so the flux 
of guided CaF molecules should be about IQ^^ per abla- 
tion shot. We see that, in combination with a buffer gas 
source, the AG guide can deliver very high numbers of 
molecules in high-field seeking states and is an attractive 
complement to the methods already developed for ex- 
tracting low-field seeking molecules from the buffer gas 
cell [Ml El]. 

It is interesting to consider whether the DC mode of 
guiding investigated here might also be useful for trans- 
porting high-field seekers from a buffer gas source. As no 
switching is required, such a guide is far easier to imple- 
ment than the AG case. To test whether the molecules 
can be guided for long periods of time, we simulated the 
molecular trajectories for times up to 500 ms. We found 
that, given enough time, all the molecules observed in 
the current experiment will crash into the electrodes even 
though some complete a large number of oscillations in 
the guide. For example, if the molecules spend 9 ms in 
the guide, instead of the current 1.7 ms, the number of 
guided molecules falls by a factor of 10, and if they spend 
50 ms in the guide the number falls by a factor of 50. We 
see that the molecules are not guided indefinitely by the 
static field but are gradually lost from the guide. This 
is in contrast to the AG guide whose phase-space accep- 
tance does not depend on the transit time once this is 
longer than the period of the macromotion. 
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